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Abstract 

We report on the formation and development of the photonic band gap in two-dimensional 8-, 10- 
and 12-fold symmetry quasicrystalline lattices of low index contrast. Finite size structures made of 
dielectric cylindrical rods were studied and measured in the microwave region, and their properties 
compared with a conventional hexagonal crystal. Band gap characteristics were investigated by 
changing the direction of propagation of the incident beam inside the crystal. Various angles of 
incidence from 0° to 30° were used in order to investigate the isotropic nature of the band gap. 
The arbitrarily high rotational symmetry of aperiodically ordered structures could be practically 
exploited to manufacture isotropic band gap materials, which are perfectly suitable for hosting 
waveguides or cavities. 

PACS numbers: 42.70.Qs, 41.20.Jb,61.44.Br 
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I. INTRODUCTION 



Structures exhibiting photonic band gap (PBG) characteristics are useful in confining and 
guiding electromagnetic energy. Photonic crystals (PCs) are artificially engineered materials 
with spatially modulated refractive indices that are widely used for such purposes. There is 
a tremendous interest in studying their properties because of their potential application in 
design and manufacturing of new optical components and devices like wavelength division 
multiplexers, switches, light emitting diodes and lasers^-. 

Recently, structures lacking long-range translational order but with orientational order 
and higher order rotational symmetries which are not compatible with the spatial periodic- 
ity, called photonic quasi-crystals (PQCs), are also gaining attention because of their unique 
characteristics. Photonic quasicrystals have neither true periodicity nor translational sym- 
metry, however they can exhibit symmetries that are not achievable by conventional periodic 
structures. These features have recently attracted a large interest because of their poten- 
tial impact in engineering novel optical circuits. PQCs can show directive emission 5 , mode 
confinement^, superlensing? , as photonic crystals do. 

One-dimensional photonic band-gap has been observed in dielectric multilayers stacked 
according to Fibonacci series^. Dielectrics arranged according to quasi-periodic geometries 
such as octagonal (8-fold), decagonal (10-fold) and dodecagonal (12-fold) are shown to have 
two-dimensional PBG 9-12 . Penrose-tiled (10-fold rotational symmetry) PQCs were the most 
studied structures among those presented above, and there were numerous studies about 
the mechanism of formation of PBG~S, optical properties, diffraction pattern^ and multiple 
scattering^ for this geometry. An organic laser based on Penrose-tiled PQCs was also 
demonstrated^. Formation of complete PBG in 12-fold symmetric dodecagonal PQC was 
numerically and experimentally studied by Zoorob et al.— . 

The lack of periodicity renders the study of PQCs very complex and computationally 
demanding. Although some concepts developed for periodic PCs can be used for the analysis 
of the properties of PQCs, a rigorous extension of a Bloch-type theorem (and associated tools 
and concepts) does not exist. Previous methods include a supercell approach in a plane- wave 
expansion technique to approximate the response of the infinite aperiodic lattice 9 - 1 ^, or the 
use of Archimedean-like tilings with properties similar to those of photonic quasi-crystals^. 
Recently, a method was proposed to computing the spectra and the eigenstates of a PQC by 
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directly solving Maxwell equations in a periodic unit cell of a higher-dimensional lattice^. 

In the geometries mentioned above, the tiling is previously calculated by matching rules 
or inflation/deflation algorithms or using projection methods from a hypercubic lattice^. 
Well-known tilings are the octagonal Ammann-BeenkerJ^, the decagonal Penrose^, and 
the dodecagonal pattern based on the Stampfli rule^, which represent the quasiperiodic 
structures object of the present study. 

The interference patterns formed by multiple-beam interferometry provide another way to 
obtain highly symmetric quasicrystalline patterns^. Recently, holographic techniques based 
on interferential methods are used to fabricate PQCs with very high rotational symmetry, 
the fold being determined by the number of interfering beams. The PBG properties of 12- 
fold symmetric quasi-crystal patterns formed by double-beam multiple exposure holography 
were studied theoretically by Gauthier et al.— . 

The method of single-beam computer-generated holography has been also successfully 
used to fabricate PQC structures with up to 23-fold rotational symmetry^. These techniques 
can be used to produce complicated two-dimensional geometries with ease. Very recently, 
a three-dimensional Penrose-type PQC fabricated through this method was also reported^. 
In considering holographic lithography, it is important to have geometries possessing op- 
timum band-gap properties at low refractive index contrast. A conceivable application is 
the combination of this versatile technology and soft materials like Polymer Dispersed Liq- 
uid Crystals for the realization of large area, high quality, low cost optical devices with 
switchable properties^. 

In two-dimensional structures, a full PBG for both transverse magnetic (TM, E field 
orthogonal to the crystal plane) and transverse electric (TE, H field orthogonal to the crystal 
plane) polarization is possible for high values of the index contrast only. Nevertheless, the 
minimum refractive index contrast AnM- at which a partial band gap starts to appear, for a 
specifically polarized electromagnetic mode, is an important parameter for the development 
of photonic devices such as waveguides. This is very much dependent on the particular 
geometry used. 

The minimum value for opening a partial band gap in a periodic (triangular) structure was 
calculated using a plane wave technique by Matthews et al.— . For the optimal configuration 
of dielectric circular rods in air (with a filling factor - the ratio between total and rods 
occupied area - r\ = 0.3), the critical value of the dielectric permittivity is 1.73 for the TM 
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FIG. 1. The periodic and aperiodic geometries under study: (a) hexagonal, (b) dodecagonal, 
(c) octagonal, and (d) decagonal (Penrose), (e): schematic diagram of the experimental set-up 
(not in scale). Horn antennas are used to transmit and receive the microwave radiation. Data 
are collected using a vectorial network analyzer computer-controlled. The size of each structure 
(periodic or aperiodic) is 40 x 14 x 60cm 3 

band gap (corresponding to An = 0.31), much larger for the TE band gap. This study 
should rule out any attempt to develop PBG devices based on periodic structures with a 
lower index contrast. 

In the case of quasi-crystal geometries, a novel method based on density waves was used to 
determine the polarized band gap for n — fold rotational symmetrical structures as a function 
of the index contrast 29 . This method works well for TM polarization only, for which the 
density function for the optimal configurations tends to have smoother features (on the scale 
of the wavelength) compared to the case of TE polarization. At low contrast, quasicrystalline 
structures of high symmetry tend to have larger band gaps than the crystalline ones, but 
smaller band gaps at high contrast. The sixfold (hexagonal) crystalline structure yields 
the largest gap for any value of the index contrast approximatively higher than 1. The 
critical value for the optimized quasi-periodic photonic crystals is estimated to be close to 
An = 0.22. 

Very recently, a detailed numerical study on the band-gap formation at low index- 
contrasts for both TE and TM polarization was carried out by Zito et al.— . Using Finite 
Difference Time Domain (FDTD) simulations, authors showed that, more than the degree 
of rotational symmetry, the difference in the tiling geometry might dramatically affect the 
existence and behavior of the band-gap. 
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Another critical parameter that comes into play for the realization of devices is the 
isotropy of the PBG. For some applications, such as light emitting diodes, this property may 
be desirable, even if the size of the full gap is slightly reduced. The periodic structures with 
square or triangular lattices have anisotropic band-gap properties because of the anisotropy 
of the Brillouin Zone (BZ). Significantly, optimized quasicrystal gaps are more isotropic than 
those of crystals, for all contrasts, due to their disallowed rotational symmetries. This is due 
to the fact that their effective BZs are more circular than the BZs of the periodic structures, 
which translates in reduced frequency undulations (that is, band gap variations over the 
wavevector in the BZ). 

Experimental investigations in this direction were performed by Bayindir et al.— in the 
microwave regime and by Hase et al.— in the far infrared region, based on octagonal and Pen- 
rose quasicrystals. Both studies reported the appearance of a PBG having almost isotropic 
properties in aperiodic lattices consisting of dielectric rods in air and for electric field parallel 
to the rods. 

We present here a detailed and systematic analysis of the band-gap isotropy of photonic 
crystal and quasi-crystal structures having low refractive index contrast and for both TM 
and TE field polarization. We consider four different geometries for comparison, the periodic 
hexagonal pattern with 6 — fold symmetry (Fig. Ufa)), and quasi-crystalline geometries with 
12 — fold symmetry (dodecagonal), 8 — fold symmetry (octagonal) and 10 — fold symmetry 
(Penrose decagonal) as shown in Fig. QJb), Fig. [Tfc) and Fig. QJd) respectively. The PBG 
properties are studied both numerically and experimentally for two different index contrasts 
close to the critical values. 

II. COMPUTATIONAL AND EXPERIMENTAL METHODS 

The photonic crystals studied consist of infinitely long dielectric cylindrical rods in air 
placed on the vertices of tiles in the corresponding geometry. The filling factor rj is set to 
be the same for all the structures under study, and is equal to 0.23. The geometries used 
are shown in Fig. [TJ They are designed to have approximately 400 rods in an area of ~ 40a 
x 14a, where a is a characteristic length of the same order of the tile side length of each 
structure. 

As discussed in the introduction, from the computational point of view each quasi- 
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crystalline structure presents a challenge in obtaining the information on the photonic band- 
gap since the lack of translational symmetry prevents the rigorous use of the Bloch theorem 
in the calculations. This difficulty has hampered the use of supercell or similar techniques 
in the modeling of real aperiodic structure. 

The Finite Difference Time Domain (FDTD) technique is useful in this respect, since 
it can provide an alternative, fast and accurate method to simulate in the real space the 
propagation of the electromagnetic waves through a finite portion of a quasi-crystal without 
recurring to any approximant. 2D FDTD method with uniaxial perfectly matched layer 
(PML) boundary conditions (along x- and y-directions) is employed to obtain transmission 
characteristics as a function of frequency of the incident radiation, propagation direction and 
polarization. In each simulation, the source of excitation consists of a time-pulsed Gaussian 
beam (with 250a < b < 800a, where b is the confocal parameter) placed outside the crystal 
structure and impinging on it with different incidence angles. The field components after 
propagation through the crystal are collected using a detector (time monitor) placed on 
the other side of the crystal. The Fourier transform of these data gives the transmission 
properties as a function of frequency. Both TM and TE polarizations were analyzed in the 
simulations. 

The experiment was designed and carried out in the microwave regime (9 — 20GHz). 
Cross-linked polystyrene (Rexolite) and polytetrafluoroethylene (Teflon) having radius 
0.64cm and length 60cm were used to fabricate the cylindrical rods. These materials 
show a dielectric constant of 2.56 and 2.1 respectively in the frequency region of interest 
and a relatively low dissipation. Loss tangent values are in the range between 10~ 3 and 
10 . In order to build up the structure, circular holes with the designed geometries have 
been drilled onto two support dielectric plates, fixed 60cm apart and then filled with the 
rods. The characteristic length a is chosen to be of the order of 1cm, so that the first PBG 
appears in the region of 10 — YlGHz for all structures under test. 

Two high-gain horn antennas acting as transmitter and receiver and connected to a two- 
port vectorial network analyzer (VNA) HP8720C have been used to obtain the transmission 
characteristics of the crystals. To change the field polarization, both horn antennas are 
rotated by 90°. The distance between transmitter and receiver is set to 3m, so that the 
polarized signal can be considered a plane wave. In the frequency region of interest, the 
signal wavelength is much smaller than the crystal height in the z-direction, therefore one 
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FIG. 2. Experimental and calculated transmission spectra of the hexagonal [(a) and (b)], dodecago- 
nal [(c) and (d)], octagonal [(e) and (f)], and decagonal [(g) and (h)] photonic crystals respectively, 
for TM polarization. Curves of different colours correspond to different angles of incidence as 
indicated in the graph. 

can safely assume that the response has a two dimensional character in the crystal (x — y) 
plane. Before each measurement the VNA is calibrated in absence of the material under 
study. Transmission curves are then obtained by introducing the sample in between the two 
horn antennas. 

In order to study the isotropic nature of the PBG, the transmission spectra are obtained 
as a function of incidence angle for all geometries under study. To record the transmission 
as a function of direction of propagation of the incident radiation, the crystal is rotated 
in 5° steps in respect to the normal direction while keeping the position of the antennas 
unchanged. Because of the n — fold rotational symmetry of the structures under test, the 
crystal/ quasi- crystal properties need only be examined over the 5° — 180°/n range since they 
repeat themselves for any propagation angle outside these degree values. For the sake of 
clarity, not all angles will be shown in the spectra of the different structures. Measurements 
were also carried out as a function of crystal thickness, for normal transmission only. 



III. RESULTS AND DISCUSSION 



The transmission characteristics are measured as a function of incidence angle and com- 
pared with simulation results for all geometries under study, using the two materials with 
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different index contrast (An = 0.6 and 0.44 for Rexolite and Tenon respectively) and for 
both TM and TE polarization. 

A. Refractive index contrast 0.60 
TM polarization: 

Fig. [2] shows the experimental and calculated transmission spectra for all the structures un- 
der study with an index contrast of 0.60 for various angles of propagation. The experimental 
results for the hexagonal photonic crystal are presented in Fig. |2fa) whereas Fig. [2(b) shows 
the corresponding simulation data. In all curves there is a strong attenuation, from around 
30 to AOdB and more, in the transmission of the electromagnetic waves through the crystal 
in the region where the first PBG is present. From the graphs it is clearly observed that 
the transmission spectra change drastically as the angle of propagation is changed from 0° 
to 30°. In the case of normal signal incidence (0°), the region of low transmission (photonic 
band-gap) is centered at about 11GHz and spans for ~ 2.6GHz, with almost no change for 
an angle 10°. However, as the angle is increased to 20°, changes are clearly visible. The cen- 
ter of the PBG is shifted to YlGHz whereas its width becomes larger (~ 3.3GHz). For an 
angle of 30°, the widening of the PBG is even stronger. Attenuation in the PBG frequency 
region increases as the incidence angle does, because of the stronger diffraction light under- 
goes due to the longer path inside the finite crystal. For the first bandgap, experimental 
results matches very well the numerical data reported in Fig. [2(b). At higher frequencies, 
a second bandgap appears with less pronounced characteristics but with similar angular 
dependence. Also in this case there is a fairly good agreement between measurements and 
simulations. 

The transmission characteristics of a dodecagonal PQC structure are shown in Figs. 
|2jc) and (d). A well-pronounced PBG, characterized by a signal attenuation larger than 
AOdB and extending approximately from 11 to 12.QGHz, is observed above the valence 
band (lower band) frequencies. At higher frequencies (upper bands) the transmitted power 
decreases, possibly with the presence of smaller and shallower bandgaps up to the maximum 
measurement frequency of 20GHz. The experimental (Fig. E(c)) and simulation (Fig. Efd)) 
angular dependence puts in evidence that in this case the PBG between the valence and 
conduction bands is nearly isotropic. Transmission spectra are only slightly affected by the 
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change of the angle of propagation from 0° to 15°. 

Similar features are observed for the other aperiodic photonic crystals. In the case of 
the 8 — fold and 10 — fold PQCs, the distribution of cylinders has a mirror symmetry with 
respect to the line of 22.5° and 18° respectively, and the data are reported in the angular 
range - 20° 

The numerical and experimental results for the case of the octagonal geometry are pre- 
sented in Figs. 12(e) and [2(f) respectively. The PBG is centered at 11.5GHz for angles from 
0° to 10° whereas it is shifted to ~ 11. 8G Hz for incidence at larger angles. There is also a 
small variation in the width of the PBG, from lAGHz to ~ 2GHz. 

Data obtained for the Penrose tiled quasicrystal are shown in Figs. 12(g) and (h). From 
both measurements (Fig. (2(g)) and simulations (Fig. (2(h)) one can clearly see that the 
bandgap is not as deep and wide as in the other cases, nevertheless its position and width 
remains almost the same for all angles. 

As in the periodic case, the transmitted signal shows an increasing attenuation in the PBG 
as a function of the incidence angle for all the aperiodic structures. Not surprisingly, the 
geometry with 12 — fold symmetry is the less sensitive to the change in the signal propagation 
direction. This is likely related to the fact that the bandgap for the dodecagonal pattern 
may be associated with a more short range ordering of the dielectric scattering centers in 
comparison with the other aperiodic structures^. 

In order to better determine the directional variation of PBG, we measured the upper 
and lower PBG boundaries as a function of the propagation angle. The band gap edges 
are defined as the frequency values for which attenuation reaches 15dB. Data in Fig. [3] 
summarize the results found from the observation of the angular dependence. The hexagonal 
geometry (Fig. E(a)) clearly shows a strong dependence of the upper and lower frequency 
edges as a function of angle, whereas the PBG width seems to be less affected. The results 
for the dodecagonal structure, instead, clearly indicates (see Fig. 13(b)) that its response 
is highly isotropic, with very small variations of the PBG width and edges at different 
angles from 0° to 15°. The properties for the octagonal geometry lie somehow in between, 
since it shows less isotropy compared to the dodecagonal geometry and a noticeable width 
dependence in respect with the hexagonal displayed in Fig. [3(c) between 0° and 

20°. Penrose geometry also seems to have quite isotropic, but narrower PBG (Fig. [3(d)) 
in the same incidence angle range. If the transmission at normal incidence is considered, 
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FIG. 3. Variation of PBG for TM polarization as a function of angle for the different geometries: 
(a) hexagonal, (b) dodecagonal, (c) octagonal and (d) decagonal. In each graph, the black and the 
red curves indicate the lower and upper frequency edges of the band gap respectively. 

the hexagonal crystal shows the widest PBG compared to all other structures under study. 
The PBG for the periodic crystal is around 2.6GHz. In the case of aperiodic structures, 
this value decreases from 2.2GHz (octagonal) to 2GHz (dodecagonal) down to 1.1GHz 
(Penrose). 



TE polarization: 

The results obtained for the hexagonal and dodecagonal geometries are compared in Fig. 
H] for normal incidence only. The hexagonal geometry (Fig. IHa)) clearly shows a PBG 
~ 1.7GHz wide centered at about 11.1GHz. The experimental (black curve) and calculated 
(red curve) results are in good agreement. The dodecagonal structure (Fig. 0Jb)) displays 
a very shallow bandgap measured approximately in the region 12 — ISGHz (black curve), 
however only partially reproduced by numerical data (red curve). 

Simulation of the response to a TE polarized wave shows instead that the transmission 
characteristics for the Penrose and octagonal structures are nearly featureless^ and are not 
reported here. 



B. Refractive index contrast 0.44 



The results presented in this case are for the hexagonal and dodecagonal cases only. The 
case of 8 — fold and 10 — fold were numerically examined in a previous paper—, where we 

10 



E 



(a) 


- 


(b) 


^ Experiment 






Simulation 







-JU H 1 , 1 , 1 , 1 , 1 , 1 H 1 , 1 , 1 , 1 , 1 , 

10 12 14 16 18 20 10 12 14 16 18 20 
Frequency (GHz) Frequency (GHz) 

FIG. 4. The transmission curves for (a) hexagonal and (b) dodecagonal geometries for normal 
incidence and TE polarization. Both experimental and simulated results are shown in each graph 
(red and black curves respectively). 
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FIG. 5. Experimental and calculated transmission spectra of the hexagonal [(a) and (b)], and 
dodecagonal [(c) and (d)] photonic crystals respectively, with an index contrast of 0.44 for TM 
polarization. Curves of different colours correspond to different angles of incidence as indicated in 
the graph. 

found that these geometries does not show any clear bandgap behavior in the frequency 
region of interest for index contrast values below 0.6 and for both polarizations. 

For TM polarization, the observed features are similar to the case of the higher refractive 
index contrast presented above. In this case too, the hexagonal crystal shows a clear variation 
in the PBG properties as the angle of propagation is varied, as seen in Figs. |5]^a) (experiment) 
and^b) (simulation). Similar results for the dodecagonal crystal are presented in Figs. El^c) 
and^d) respectively. The PBG frequency shift is reduced in comparison with the hexagonal 
crystal but the PBG is much narrower. A graph showing the variation of the upper and 
lower band-edge frequencies as a function of angle is also plotted in this case and shown 
in Fig. |6] for both the 6 — fold and 12 — fold symmetry (in the range of — 30° and 
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FIG. 6. Variation of PBG as a function of angle for (a) hexagonal and (b) dodecagonal geometries, 
for TM polarization. In each graph, the black and the red curves indicate the lower and upper 
frequency edges of the band gap respectively. 
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FIG. 7. Experimental transmission curves as a function of sample thickness for (a) hexagonal, 
(b) dodecagonal, (c) octagonal, and (d) decagonal geometries with index contrast 0.60 and TM 
polarization. 

— 15° respectively). As expected due to the lower index contrast, using Teflon cylinders 
the observed PBGs are never as pronounced as in the case of Rexolite. 

For TE polarization, the PBG features shown in the spectra were so weak to render the 
study of its angular dependence meaningless. Once again, this is in agreement with the 
previous simulation work 30 . 



C. Effects of thickness 

To understand how sensitive is the PBG to the size of a photonic quasi-crystal with 
a low index contrast, transmission spectra were also measured by removing the cylinder 
rows simultaneously from both sides of the structure in the signal propagation direction. 
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Results are shown in Fig. [7] for all aperiodic geometries with an index contrasts of 0.6 for 
normal incidence and in the bandgap region only. The hexagonal case is also reported for 
comparison. Graphs clearly show that the dodecagonal pattern is the less affected by the 
reduction in size, presenting an appreciable bandgap (15dB attenuation) even when the 
slab thickness decreases from 14a to 8a. Similar results, but obviously less pronounced, are 
also obtained for Teflon (lOdB attenuation). This is a further confirmation that the 12-fold 
geometry is the most robust amongst the different structures under investigation because 
of the stronger role played by the short-range interactions, and it is therefore preferable for 
the development of very compact photonic devices. 



IV. CONCLUSIONS 

It is well known^i that for transverse magnetic modes bandgaps can be easily achieved in 
photonic crystals made of isolated dielectric rods, whereas connected lattices favour trans- 
verse electric gaps. The same happens in photonic quasi-crystals, where we have experi- 
mentally shown that in the dielectric-in-air configuration two-dimensional TM bandgaps are 
clearly visible whereas the response of the aperiodic structures to a TE polarised signal is 
very weak or almost featureless. Besides that, results put in evidence that TM bandgaps in 
PQC are possible even with a very low-index-contrast, in agreement with numerical stud- 
ies carried out in Ref. Q. In particular, well-pronounced bandgaps are present for the 
dodecagonal geometry, as theoretically expected because of the higher symmetry. 

Also, we did observe that gaps in quasicrystalline geometries are more isotropic - although 
narrower - than those in periodic crystals, due to their disallowed, non-crystallographic, 
rotational symmetries. That is, the position and width of the PBG are almost independent 
of the incident angle of the light, contrarily to their periodic counterparts, where gaps of 
different directions may appear at different frequencies because of non-spherical first BZ. 
For specific applications, like light emitting diodes or waveguides, very isotropic PBGs may 
be desirable, even if the size of the full gap is slightly reduced. 

Another important feature of PQCs is that the existence of the gaps is governed by 
the short-range environment. This is particularly evident in the case of the dodecagonal 
geometry, where the PBG is robust even for a significant reduction in size in the propagation 
direction. 
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In conclusion, this study confirms that quasicrystalline structures having long-range ori- 
entational order forbidden for periodic systems are promising candidates as PBG materials. 
Photonic crystals based on aperiodic specific geometries present extremely interesting fea- 
tures that cannot be achieved in the periodic case. The low index contrast allows the use 
of versatile and low-cost technologies like holographic lithography combined with soft ma- 
terials for the development of compact devices with switchable properties for an all-optical 
ultrasmall integrated circuitry. 
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